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Abstract 

Extending results of Oh-Zumbrun and Johnson-Zumbrun for parabolic conservation 
laws, we show that spectral stability implies nonlinear stability for spatially periodic 
viscous roll wave solutions of the one-dimensional St. Venant equations for shallow 
water flow down an inclined ramp. The main new issues to be overcome are incomplete 
parabolicity and the nonconservative form of the equations, which leads to undiffer- 
entiated quadratic source terms that cannot be handled using the estimates of the 
conservative case. The first is resolved by treating the equations in the more favorable 
Lagrangian coordinates, for which one can obtain large-amplitude nonlinear damping es- 
timates similar to those carried out by Mascia-Zumbrun in the related shock wave case, 
assuming only symmetrizability of the hyperbolic part. The second is resolved by the 
observation that, similarly as in the relaxation and detonation cases, sources occurring 
in nonconservative components experience greater than expected decay, comparable to 
that experienced by a differentiated source. 

1 Introduction 

Roll waves are a well-known hydrodynamic instability occurring in shallow water flow down 
an inclined ramp, generated by competition between gravitational force and friction along 
the bottom. These can be modeled as periodic traveling- wave solutions of the St. Venant 
equations for shallow water flow, which take the form of hyperbolic or parabolic balance 
laws; see [Dj INH IN2| for detailed discussions of existence in the inviscid and viscous case. 
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The spectral and linear stability of roll waves has been studied for the inviscid St. Venant 
equations in [Nlj and the viscous St. Venant equations in |N2| . However, up to now, the 
relation between spectral, linearized, and nonlinear stability has remained an outstanding 
open question. In this paper, extending recent results of [OZ4t I JZ3|, IJZ4] in the related 
conservation law case, we settle this question by showing that spectral implies linearized 
and nonlinear stability. 

This opens the way to rigorous numerical and analytical exploration of stability of roll 
waves and related phenomena via the associated eigenvalue ODE, a standard and numer- 
ically and analytically well-conditioned problem. At the same time, it gives a particularly 
interesting application of the techniques of |OZ4t [JZ3l IJZ4| . For, roll waves, by numerical 
and experimental observation, appear likely to be stable, at least in some regimes. In the 
conservation law case, by contrast, periodic waves so far appear typically to be unstable 



1.1 Equations and assumptions 

Consider the one-dimensional St. Venant equations approximating shallow water flow on 
an inclined ramp: 



(1.1) 



h t + (hu) x = 0, 
(hu) t + (h 2 /2F + hu 2 ) x = h-u 2 + u(hu x ) x , 



where h represents height of the fluid, u the velocity average with respect to height, F is the 
Froude number, which here is the square of the ratio between speed of the fluid and speed of 
gravity waves, v = Re -1 is a nondimensional viscosity equal to the inverse of the Reynolds 
number, the term u 2 models turbulent friction along the bottom, and the coordinate x 
measures longitudinal distance along the ramp. 
In Lagrangian coordinates, these appear as 

Tt — u x = 0, 

(12) 

u t + {{2F)- 1 t- 2 ) x = 1-tu 2 + v(t- 2 u x ) x , 

where r := h^ 1 and x now denotes a Lagrangian marker rather than physical location. We 
will work with this form of the equations, as it is more convenient for our analysis in several 
ways. (Indeed, for the large-amplitude damping estimates of Section I4.1|, it appears to be 
essential in order to obtain quantitative bounds on amplitude; see Remark 14.31 ) 
Denoting U := (r, it), consider a spatially periodic traveling-wave solution 

(1.3) U = U(x-ct), 

of (II. 2b of period X and wavespeed c satisfying the traveling-wave ODE 

— cr — u = U, 

' ' ' ' ' -cv! + ((2F)- 1 r- 2 ) / = 1 - tu 2 + v(t- 2 U ')', 
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Integrating the first equation of (jl.4p and solving for u = u(t) := q — cr, where q is 
the resulting constant of integration, we obtain a second-order scalar profile equation in r 
alone: 

(1.5) cV + ((2F)-V- 2 )' = l-r(q- ct) 2 - cv{t~ 2 t') '. 

Note that nontrivial periodic solutions of speed c = do not exist in Lagrangian coordinates, 
as this would imply u = q, and fjl .5[) would reduce to a scalar first-order equation 

(1.6) r' = Fr\rq 2 -l), 

which since it is scalar first-order has no nontrivial periodic solutions, even degenerate 
ones (e.g., homoclinic or heteroclinic cycles) that might arise in the singular c — > limit. 
Rather, there appears to be a Hopf bifurcation as c approaches some minimum speed for 
which periodics exist; see [N2| . Section 4.1 and Fig. 1, Section 4.2.3. 

It follows then that periodic solutions of (|l,5p correspond to values (X,c,q,b) £ M 5 , 
where X, c, and q denote period, speed, and constant of integration, and b = (61,62) 
denotes the values of (r, r') at x = 0, such that the values of (r, t') at x = X of the solution 
of (jl.5p are equal to the initial values (61, 62)- 

Following [Sell IOZ31 IQZil lJZ3l iJZl] , we assume: 

(HI) f > 0, so that all terms in ([L2} are C K+1 , K > 3. 

(H2) The map H : M. 5 ->• M 2 taking (X,c,q,b) ^ ( r , r') (X, c, 6; X) - b is full rank at 
(X,c, b), where (r, t')(-;-) is the solution operator of f)1.5|) . 

By the Implicit Function Theorem, conditions (H1)-(H2) imply that the set of periodic 
solutions in the vicinity of U form a smooth 3-dimensional manifold 

(1.7) {UP(x -a- c{/3)t)}, with aelJeM 2 . 

Remark 1.1. The transversality condition (H2) could be replaced by the more general as- 
sumption that the set of periodic solutions in the vicinity of U form a smooth 3-dimensional 
manifold ()1.7|) . However, it is readily seen in this context that (H2) is then implied by the 
spectral stability condition (D3) of Section [l.l.2j that is, transversality is necessary for our 
notion of spectral, or Evans, stability. This situation is reminsiscent of that of the viscous 
shock case; see, for example, |ZH[ S 1.2.3], or [MaZ3, Zl . 

Remark 1.2. Note that (jl-2j) is of 2 x 2 viscous relaxation type 



(1.8) V, -^D(U)U r ), : ( g( ^) '/(f) !-r<r 



where q u = —2ut < for solutions u > progressing down the ramp. T 

u 3 
2 



solutions are stable so long as the subcharacteristic condition 



< 



rus, constant 
is satisfied, 



or F < 4. When the subcharacteristic condition is violated, roll waves appear through 
Hopf bifurcation as parameters are varied through the minimum speed c mm = , 1 ; see 



/Ft 3 ' 



Appendix [Cj For v = 0, violation of the subcharacteristic condition is associated with 
subshocks and the appearance of discontinuous roll waves observed by Dressier [D]; see [JK| 
for related, more general, discussion. 
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Remark 1.3. The limit v — >• represents an interesting singular perturbation problem in 
which the structure of the profile equations simplifies, decoupling into fast and slow scalar 
components, and converging to inviscid Dressier waves (Dj IN1| in an appropriate regime 
[N2j . This would be an interesting setting in which to investigate the associated spectral 
stability problem. Another interesting limit is Hopf bifurcation from the constant solution 
occurring at minimum speed of existence |.\2j , treated here in Section [Cj see Remark IC. II 

1.1.1 Linearized equations 

Making the change of variables x — > x — ct to co-moving coordinates, we convert (|1.2|) to 

TV — ct x — u x = 0, 

(19) 

u t - cu x + ((2F)- 1 r- 2 ) a . = 1 - tu 2 + v(t~ 2 u x ) x , 

and the traveling-wave solution to a stationary solution U = U(x) convenient for stability 
analyis. 

Writing (11. 9ft in abstract form 

(1.10) U t + f(U) x = {B(U)U x ) x + g{U) 
and linearizing (|1 .9j) about £/(•), we obtain 

(1.11) v t = Lv:= (d x Bd x -d x A + C)v, 
where the coefficients 

A := df(U) - ( d B(U)W* = (_,_3 (F -7_ 2v - Ux) ~_l 

(1.12 N 



*==*<*>=(!! JU). ^*m=(^_L 

are periodic functions of x. As the underlying solution U depends on x only, equation (jl.lip 
is clearly autonomous in time. By separation of variables, therefore, decomposing solutions 
into the sum of solutions of form v(x, t) = e xt v(x), where v satisfies the eigenvalue equation 
(L — X)v = 0, or, equivalently, by taking the Laplace transform, we may reduce the study 
of stability of U to the study of the spectral properties of the linearized operator L. 

As the coefficients of L are X-periodic, Floquet theory implies that its spectrum is 
purely continuous. Moreover, its spectral properties may be conveniently analyzed by Bloch 
decomposition, an analog for periodic-coefficient operators of the Fourier decomposition of 
a constant-coefficient operator, as we now describe. 

1.1.2 Bloch decomposition and stability conditions 

Following [Gl EH IS2l [S3] , we define the family of operators 

(1.13) Lt. = e-^Le^ = (d x + it)B{d x + i£) - (d x + %)A + C 
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operating on the class of I? periodic functions on [0, X\\ the (L 2 ) spectrum of L is equal to 
the union of the spectra of all with £ real with associated eigenfunctions 

(1.14) w(x,Z,\):=e* x q(x,Z,\), 

where q, periodic, is an eigenfunction of L^. By standard considerations [N2] F1 the spectra 
of consist of the union of countably many continuous surfaces Aj(£). 

Without loss of generality taking X = 1, recall now the Block representation 

(1.15) u(x)= (^) £ e**ufox)d£ 

of an L 2 function u, where u(£, x) := ^ fc e 27rifca: n(£ + 2-7r/c) are periodic functions of period 
X = 1, u(-) denoting with slight abuse of notation the Fourier transform of u in x. By 
Parseval's identity, the Bloch transform u(x) — > u(£,x) is an isometry in L 2 : 

(1-16) IMIl 2 (x) = \\u\\l 2 (c,l 2 (x)), 

where L 2 (x) is taken on [0, 1] and L 2 (^) on [— tt, it]. Moreover, it diagonalizes the periodic- 
coefficient operator L, yielding the inverse Bloch transform representation 

(1-17) e Lt u = (i-) £ e^e^u (£,xR 

relating behavior of the linearized system to that of the diagonal operators L^. 

Following |JZ4| . we assume along with (H1)-(H2) the strong spectral stability conditions: 
(Dl) <t{L{) C {ReA < 0} for £ ^ 0. 

(D2) Rea{L/:) < -6\£\ 2 , 9 > 0, for £ E R and |£| sufficiently small. 
(D3') A = is an eigenvalue of Lq of multiplicity 20 

As shown in |N2j . (H1)-(H2) and (D1)-(D3') imply that there exist 2 smooth eigenvalues 

(1.18) A i (0 = -m^ + o(|e|) 

of bifurcating from A = at £ = 0; see Lemma 12 . 1 1 b elow . 

Loosely following |JZ4| . we make the further nondegeneracy hypotheses: 

(H3) The coefficients aj in (|1.18p are distinct. 

(H4) The eigenvalue of Lq is nonsemisimple, i.e., dimker Lq = 1. 
The coefficients aj may be seen to be the characteristics of an associated Whitham averaged 
system 

fllQl M((3) t + G((3) x = 0, 

n(p) t + (c(p)n(p)) x = o 



For example, the characterization [G] of spectra as the zero set of an associated Evans function. 
2 The zero eigenspace of Lo, corresponding to variations along the 3-dimensional manifold of periodic 
solutions in directions for which period does not change Scl ( JZ4 , is at least 2-dimensional by linearized 
existence theory and (H2). 
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linearized about the values of M, G, c, associated with the background wave u, where M 
is the mean of r over one period and F the mean in the r-coordinate of a certain associated 
flux, c is wave-speed, and VL frequency of nearby periodic solutions, indexed as in (jl.7p by 
/9 £ E 2 ; see [N2l IOZ31 lOZi] System fpT~T9^ formally governs slowly modulated solutions 

(1.20) u(x,t) =ul 3{£X > £t \y(x,t)) +0(e), e^O 

presumed to describe large spatio-temporal behavior x, t 3> 1, where Cr(-) as in (jl.7p 
parametrizes the set of nearby periodic solutions, Q = ty x , and c = — ^t/^x- 

Thus, (Dl) implies weak hyperbolicity of the Whitham averaged system fjl . 191) (reality 
of a,j), while (H3) corresponds to strict hyperbolicity. Condition (H4) holds generically, and 
corresponds to the assumption that speed c is nonstationary along the manifold of nearby 
stationary solutions; see Lemma 12. ill Condition (D2) corresponds to "diffusivity" of the 
large-time (~ small frequency) behavior of the linearized system, and holds generically given 
(H1)-(H4), (Dl), and (D3')d Condition (D3') also holds generically, and can be verified by 
an Evans function computation as described in [Nl| . As discussed in [OZl] ISeTj IJZ3} |JZ4] . 
conditions (D1)-(D3') are conservation law analogs of the spectral assumptions introduced 
by Schneider in the reaction-diffusion case |SH |S2} [S3] . 



1.2 Main result 

Theorem 1.1. Assuming (H1)-(H4) and (D1)-(D3 J ), let U = (f,u) be a traveling-wave 
solution (jl.3p of (|1 .2j) satisfying the derivative condition 

(1.21) uu x < F~ l . 

Then, for some C > and tp G W K,oc (x,t), where K > 3 is as in (HI) 

\\u -u(--i>- ct)\\ LP (t) < c(i + ty^-y^wu - u\\ L , nH K\ t=0 , 

(1.22) \\U -U(--iP- ct)\\ HK {t) < C(l + t)-\\\U - C7|| iln ^k=o ; 

\\WM\\w«+^ <C(l + *)~^ (1 ~ 1/p) ll^-^llLin^lt=o, 

and 

(1.23) \\U - U(- - ct)\\ L ~(t), < C\\U - U\\ L i nH K\t=o 

for all t > 0, p > 2, for solutions U of (II. 2b with \\U — U\\ L i nH ic\t=o sufficiently small. In 
particular, U is nonlinearly bounded L 1 n H — > L°° stable. 

Theorem 11.11 asserts not only bounded L 1 n H K — > L°° stability, a very weak notion of 
stability, but also asymptotic convergence of U to the modulated wave U(x — ip(x,t)). 

3 Here, we follow the formalism and notation of JOZ3I [OZ4] , 

4 The case that (H4) is violated may be treated as in [JZ3| . 

5 This amounts to nonvanishing of bj in the Taylor series expansion Aj(£) = —iaj£ — bj£ 2 guaranteed by 
Lemma [2T1 given (H1)-(H4), (Dl), and (D3'). 
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Remark 1.4. With further effort, it may be shown that the results of Theorem II . 1 1 extend 
to all 1 < p < oo using the pointwise techniques of |OZ2| ; see discussion, [JZ3, IJZ4] . 

Remark 1.5. The derivative condition (jl.2ip is effectively an upper bound on the ampli- 
tude of the periodic wave; see Remark 14. 2 [ As discussed in Remark 14.31 this is precisely the 
condition that the first-order part of the linearized equations (jl.lip be symmetric hyper- 
bolic (i.e., that A in (11. 12ft be symmetrizable) , and reflects a subtle competition between 
hyperbolic and parabolic effects. (The first-order part of the inviscid equations is always 
symmetric-hyperbolic, corresponding to the equations of isentropic gas dynamics with 7- 
law gas.) It is satisfied when either wave amplitude or viscosity coefficient v is sufficiently 
small. It is not clear whether this condition may be relaxed. 

We note that condition (jl.21[) is satisfied for all roll-waves computed numerically in 
[N2j . For, in Eulerian coordinates, this condition translates to h x /h < (cuF)^ 1 . Examining 
Fig. 1 of [N2j . a phase portrait in (h,h r ) for F = 6, v = 0.1, and 1.89 < c < 1.91, we see 
that all periodic orbits appear to lie beneath the line h! jh = .68, whereas (cvF)~ 1 ~ .88. 

It is straightforward using the bounds of Corollary 13.11 to show for "zero-mass" , or 
derivative, initial perturbations, that nonlinear decay rates (ll.22p - (|1.23p improve by fac- 
tor (1 + i) -1 / 2 , to the rates seen in the reaction-diffusion case |SH |JZ5] for general (un- 
differentiated) localized perturbations. In particular, the perturbed wave U then decays 
asymptotically in L°° to the background wave U with Gaussian rate (1 + i) -1 / 2 as in the 
reaction-diffusion case. Likewise, under an unlocalized initial perturbation, or, equivalently, 
the integral of a localized perturbation the difference between U and U may be expected to 
blow up at rate (1 + 1) 1//2 - this is indeed the linearized behavior- and, barring special non- 
linear structure, there seems no reason why the difference between U and the modulation 
U(- — \I/) should not blow up as well: at best it remains bounded. In the reaction-diffusion 
case, for comparison, results announced in [SSSUJ assert that U remains close to U even 
under unlocalized perturbations, and approaches the modulated wave at rate (1 + t)~ 1 / 2 in 
L°° . That is, the behavior in the conservation (balance) law case compared to that in the 
reaction-diffusion case is, roughly speaking, shifted by one derivativejj 

This reflects a fundamental difference between modulational behavior in the present, 
conservation (or balance) law setting from that of the reaction-diffusion case. Namely, 
in the reaction-diffusion case, the Whitham averaged system reduces to a single equation 
<9f($7) + d x (flc) = 0, or, equivalently, 

(1.24) y t + c(y x )V x = 0, 

where £1 := SS? X denotes frequency and c := — wave speed, and c and are related by 
the linearized dispersion relation along the family of periodic orbits (in the case considered 
by Schneider |S1| . c = 0). On the other hand, the Whitham averaged equations (|1.19p in 

6 At a purely technical level, this can be seen by the appearance of a Jordan block in the zero eigenspace 
of Lq, introducing factor f; -1 in the description of low- frequency behavior (Lemma l2.1[) . Recall that a factor 
i£ corresponds roughly to differentiation in the Bloch representation, through its relation to the Fourier 
transform. In the reaction-diffusion case, the zero eigenspace of Lq is simple, and no such factor appears. 
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the present case are a genuine 2x2 first-order hyperbolic systerd_| in *& x and wave-speed 
c, c now considered as an independent parameter; that is, they describe modulation of the 
perturbed wave in frequency ^f x and speed c, with phase shift ^ determined indirectly by 
integration of ^f x . 

Assuming heuristically (as justified at the linearized, spectral, level by the Bloch analysis 
of Section [2]) , that modulational behavior is governed by a second-order regularization of 
the first-order Whitham averaged system, we have the standard picture of behavior under 
localized perturbation as consisting of modulations in (^f x , c) given by a pair of approximate 
Gaussians propagating outward with Whitham characteristic speeds a\ and 02, hence an 
associated, much larger modulation in determined by integration in the ^f x component, 
given by a sum of approximate errorfunctions propagating with the same speeds. 

Indeed, this is exactly the description given in (|3.27[) of the principal part of the kernel 
e(x,t;y) determining ^ through (|4.23p . Likewise, the principal part of the Green function 
of the linearized equations about U is U'(x)e(x,t;y), showing that linearized behavior to 
lowest order indeed consists of a translation, or multiple of U'(x), with amplitude 




(x,t;y)(U(y,0)-U(y,0))dy; 



see the description of the Green function in Corollary 13.11 The same considerations show 
that the rate of convergence of U to the modulation U(- — \P) cannot be improved (or, in 
the case of a nonlocalized perturbation, recovered) by modulating in additional parameters 
such as wave speed c or etc. For, as indicated by the above discussion, all such modulations 
represent smaller contributions by factor (1 + i) -1 / 2 than that of the phase *S>, comparable 
rather to frequency ^f x , and thus may be ignored in consideration of blow-up vs. stability. 

This picture of modulational behavior as "filtering" by integration along a certain di- 
rection of the hyperbolic-parabolic system derived by Whitham averaging seems quite in- 
teresting at a phenomenological level, and a genuinely novel aspect of the conservation 
(balance) law case. In particular, the ^f x component direction along which the integration 
is performed is in general independent of either characteristic mode, so that the resulting 
behavior is essentially different from that exemplified by (11. 24ft of a single scalar equation 
as in the reaction-diffusion case. 

1.3 Discussion and open problems 

The extension from the parabolic conservation law to the present case involves a number 
of new technical issues associated with lack of parabolicity and nonconservative form. We 
overcome these difficulties by combining the arguments of |JZ3l IJZ4| . |N2| with those of 
|MaZ41 IZil ITZT] (real viscosity) and [MaZlj and [LRTZ1 [TZ2] (relaxation and combustion 
systems both involving nonconservative terms). 

An interesting open problem is the rigorous justification of spectral stability of roll waves 
approaching the inviscid case in the singular zero viscosity limit, extending results of |N2j. 

7 In general, the dimension of the Whitham averaged system is equal to the dimension of the manifold 
of nearby periodic solutions, modulo translations [JZ4| . 
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We hope to carry this out in future work. For related asymptotic analysis, see the study in 
[Z2J of the inviscid limit for detonations. 

Another interesting open problem is the numerical investigation of spectral stability of 
large-amplitude roll waves. In particular, it is an interesting question whether violation 
of the apparently technical "amplitude condition" (|1.2ip corresponds to actual physical 
phenomena/instability. This is not inconceivable, as (jl.2ip is needed in our argument not 
only for nonlinear iteration, but also for high-frequency linearized bounds. As the condition 
that the first-order part of the equations be symmetric hyperbolic, it may well have such 
significance- however, this is not yet clear. 

It is straightforward to extend our results to the two-dimensional small-amplitude case, 
by working in Eulerian coordinates and substituting for the present large-amplitude damp- 
ing estimate the simpler small-amplitude version of [MaZ2]: see \JZ3\ IJZ4] for the multi- 
dimensional analysis of periodic waves. However, there is some evidence that roll waves 
develop transverse instabilities in multi-dimensions [N3]. If so, this suggests the ques- 
tion whether such instability might be connected with bifurcation to multiply periodic 
waves. The extension of our stability analysis to the multiply periodic case, as suggested 
in [JZ3|. [JZ4] . would be another very interesting open problem. 



2 Spectral preparation 

We begin by a careful study of the Bloch perturbation expansion near £ = 0. 

Lemma 2.1. Assuming (H1)-(H4), (Dl), and (D3'), the eigenvalues A,(£) of are 
analytic functions and the Jordan structure of the zero eigenspace of Lq consists of a 1- 
dimensional kernel and a single Jordan chain of height 2, where the left kernel of Lq is 
spanned by the constant function f = (1,0) T , and v! spans the right eigendirection lying 
at the base of the Jordan chain. Moreover, for |£| sufficiently small, there exist right and 
left eigenf unctions qj(£, •) and qj(£, •) of L^ associated with Xj of form qj = Yl1=i Pj,k v k 
and q~j = Y^1=i Pj,kVk where {vj}j =1 and {vj} 2 j = \ are dual bases of the total eigenspace 
of L^ associated with sufficiently small eigenvalues, analytic in £, with {^(O) constant and 
v\(fd) = u'(-); £ _1 /3j,i, /?j,2 and /3j,2 ire analytic in £; and (q~j,qk) = Sj- 

Remark 2.2. Notice that the results of Lemma 12.11 are somewhat unexpected since, 
in general, eigenvalues bifurcating from a non-trivial Jordan block typically do so in a 
nonanalytic fashion, rather being expressed in a Puiseux series in fractional powers of 
The fact that analyticity prevails in our situation is a consequence of the very special 
structure of the left and right generalized null-spaces of the unperturbed operator Lq, and 
the special forms of the equations considered. 

8 This is, however, consistent with the picture of behavior as being approximately governed by a first- 
order Whitham averaged system with eigenvalue perturbation expansions agreeing to first-order with the 
associated linearized homogeneous dispersion relation [N2l IOZ31 [QZ4] . 
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Proof. Recall that Lg has spectrum consisting of isolated eigenvalues of finite multiplicity 
(HUE]. Expanding 

(2.1) L^ = Lo + i^L l -£ 2 L 2 , 
where, by (|1 . 13|) . 

(2.2) L = d x Bd x -d x A + C, L 1 = (Bd x + d x B — A), L 2 = B, 

consider the spectral perturbation problem in £ about the eigenvalue A = of Lq. 

Because is an isolated eigenvalue of Lq, the associated total right and left eigenprojec- 
tions Pq and Po perturb analytically in £, giving projection Pg and Pg [Kj. These yield in 
standard fashion (for example, by projecting appropriately chosen fixed subspaces) locally 
analytic right and left bases {vj} and {vj} of the associated total eigenspaces given by the 
range of Pg, Pg. 

Defining V = (vi, 1)%) and V = (v\, V2)*, * denoting adjoint, we may convert the infinite- 
dimensional perturbation problem (|2.ip into a 2 x 2 matrix perturbation problem 

(2.3) M 5 = M + z£Mi -£ 2 M 2 + 0(|£| 3 ), 

where := (vg , L^Vg) and (•, •) denotes the standard L 2 (x) inner product on the finite 

interval [0,X]. That is, the eigenvalues Aj(£) lying near of Lg are the eigenvalues of Mf, 
and the associated right and left eigenfunctions of L^ are 

(2.4) fj = Vwj and fj = WjV* , 

where Wj and Wj are the associated right and left eigenvectors of Mg. 

By assumption, A = is a nonsemisimple eigenvalue of Lq, so that Mq is nilpotent 
but nonzero, possessing a nontrivial associated Jordan chain. Moreover, using the fact 
that ((1,0) T ,C) = 0, where, again, (•, •) represents the L 2 (x) inner product over the finite 
domain [0, X], the function / = (1,0) T by direct computation lies in the kernel of Lq = 
(d x B*d x + A*d x + C*), we have that the two-dimensional zero eigenspace of Lq is consists 
precisely of a one-dimensional kernel and a single Jordan chain of height two. Moreover, by 
translation-invariance (differentiate in x the profile equation (II. 5p ). we have Lqu' = 0, so 
that v! lies in the right kernel of Lq. 

Now, recall assumption (H2) that H : E 5 -)• R 2 taking (X, c, q, b) (->• (r, t')(X, c, b;X)-b 
is full rank at (X,c,b), where (r, t' )(■;■) is the solution operator of f)1.5|) . The fact that 
kerLo is one-dimensional implies that the restriction H taking (b,q) 1— > u(X;b,c,q) — b for 
fixed (X,c) is also full rank, i.e., H is full rank with respect to the specific parameters 
(X, c). Applying the Implicit Function Theorem and counting dimensions, we find that 
the set of periodic solutions, i.e., the inverse image of zero under map H local to u is a 
smooth three-dimensional manifold {u^{x — a — c(/3)t)}, with a G M, /3 G M 2 . Moreover, 
two dimensions may be parametrized by (X, c) , or without loss of generality f3 = (X, c) . 
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Fixing X and varying c, we find by differentiation of (jl.5p that /* := —d s U satisfies the 
generalized eigenfunction equation 

L f* = U'. 

Thus, U' spans the eigendirection lying at the base of the Jordan chain, with the generalized 
zero-eigenfunction of Lq corresponding to variations in speed along the manifold of periodic 
solutions about U. Without loss of generality, therefore, we may take v 2 to be constant at 
£ = 0, and v x = U' at £ = 0. 

Noting as in [JZ3] the fact that, by (fl~T2|) . 

AU X = f(U) x - (d x B(U))U x = d x (f(U) x - B{U)U X ) + B{U)d x U x 

= g(U) + B(U)d x U x , 

and so by e 2 g = 0, d x e 2 = 0, we have 

{e 2 ,L 1 U') = <e 2 , (8 X B + Bd x - A)U') = {e 2 ,d x BU') = 

for e 2 := (0, 1), where (•, •) denotes L 2 (x) inner product on the interval x G [0, X], we find 
under this normalization that (|2.3h has the special structure 

<*» *=(;:)■ 

Now, rescaling (12, 3p as 

a/ € : = m-'sm^s^r 1 , 



M c = M + ^Af! + O(e 2 ), 
where M,- = M,- like the original Mj are constant and the eigenvalues rrij{^) of are 

As the eigenvalues rrij of are continuous, the eigenvalues = i^Trij are differen- 

tiable at £ = as asserted in the introduction. Moreover, by (H3), the eigenvalues Xj(0) 
of Mq are distinct, and so they perturb analytically in £, as do the associated right and 
left eigenvectors Zj and %. Undoing the rescaling (|2.7p . and recalling (|2.4p . we obtain the 
result. □ 



(2.7) 
where 

(2.8) 

we obtain 
(2.9) 
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3 Linearized stability estimates 

By standard spectral perturbation theory [K] . the total eigenprojection -P(£) onto the 
eigenspace of L^ associated with the eigenvalues Aj(£), j = 1,2 described in the previous 
section is well-defined and analytic in £ for £ sufficiently small, since these (by discreteness 
of the spectra of L^) are separated at £ = from the rest of the spectrum of Lq. By (D2), 
there exists an e > such that KAj(£) < — #|£| 2 for < |£| < 2e. With this choice of e, we 
introduce a smooth cutoff function that is identically one for |£| < e and identically 
zero for |£| > 2e, e > sufficiently small, we split the solution operator S(t) := e Lt into a 
low-frequency part 

(3.1) S^uo := (^) e^mPi^Mi, 
and the associated high-frequency part 

(3.2) S H (t)U := (i-) |" e^(l - ^(O^L/b^s)^ 

Our strategy is to treat the high- and low-frequency operators separately since, as is 
standard, the low-frequency analysis is considerably more complicated than the correspond- 
ing high-frequency analysis. That being said, we begin by deriving bounds on the solution 
operator at high frequency. 

3.1 High-frequency bounds 

By boundedness of the resolvent on compact subdomains of the resolvent set, equivalence (as 
the zero-set of an associated Evans function |N21IOZi] ) of H l and I? spectrum, Assumption 
(D2), and the high-frequency estimates of Lemma EH we have for |£| bounded away from 
zero and waves satisfying the amplitude condition (jl.2ip that the resolvent {\ — L^) v is 
uniformly bounded from H 1 — > H 1 for KA = —rj < 6 < 0, whence, by Priiss' Theorem |Prj . 
\\e L ^f\\ m <Ce- et \\f\\ H t. 

For |£| sufficiently small, on the other hand, = 1, and I — <ft(t;)P = I — P = Q, where Q 
is the eigenprojection of L^ associated with eigenvalues complementary to Aj(£), which by 
spectral separation of Aj(£) from the remaining spectra of L^, have real parts strictly less 
than zero. Applying Priiss' Theorem to the restriction of L^ to the Hilbert space given by 
the range of Q, we find, likewise, that ||e L e*(J - <t>(0)f\\m = \\e L ^Qf\\m < Ce~ et \\f\\ H i . 

Combining these observations, we have the exponential decay bound 

\\e L t\l-ct>P(0)f\\m({o,x]) < Ce- et \\f\\ mmx]) 
for 9 > as in (D2) and C > 0, from which it follows 
(3-3) \\e L ^(I - <f>P(0)d l J\\ H mo,X]) < Ce- et \\f\\ Hl+H[0tX]) 

for < / < K (K as in (HI)). Together with (|1.16p . these give immediately the following 
estimates. 
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Proposition 3.1 Q0Z4J). Under assumptions (Hl)-(H^), (D1)-(D2), and assuming the 
amplitude condition (jl.2ip holds, there exists constants 9, C > 0, such that for all all t > 0, 
2 < p < oo, < I < 2, < m < 2, we have the high-frequency estimates 

ll^ 7 (t)5i/||L 2( ,)<Ce- e '||/|| Hl+1(:E) , 
\\S U {t)d™f\\ LP{x) < Ce- Bt \\f\\ Hm+ . {x) . 

Proof. For m, I = 0, the first inequalities follow immediately by (11. 16ft and (j3.3[) . The second 
follows for p = oo by Sobolev embedding. The result for general 2 < p < oo then follows by 
IP interpolation. A similar argument applies for 1 < l,m < 2 by higher-derivative versions 
of (|3.3p . which follow in exactly the same way. □ 

3.2 Low-frequency bounds 

As noted above, analysis of the solution operator at low frequency is considerably more 
complicated than the high-frequency bounds outlined above. To aid in our analysis, we 
introduce the Green kernel 

(3.5) G I (x,t;y):=S I (t)S y {x) 
associated with S 1 , and the corresponding kernel 

(3.6) [Gl(x,t;y)] := <K0^(£)e L? %(*)] 

appearing within the Bloch representation of G 1 , where the brackets on [G^] and [5 y ] denote 
the periodic extensions of these functions onto the whole line. Then, we have the following 
descriptions of G 1 , deriving from the spectral expansion (|1.18p of near £ = 0. 

Proposition 3.2 ( [OZ4] ). Under assumptions (H1)-(H4) and (Dl)-(DS'), 

2 

[G\{x, t- y)] = m e Aj ®% (£, x)qj (£, y)* , 

(3.7) G'ix^y) = (i-) J e*<*-v\G\{x,tw)]dt 

1 r 2 

= (^) J^^^^e^qiiM^y)*^, 

where * denotes matrix adjoint, or complex conjugate transpose, qj(£, •) and q~j(£, •) are right 
and left eigenf unctions of L^ associated with eigenvalues defined in (|1.18j) . normalized 

so that (q~j, qj) = 1. 

Proof. Relation (I3.7P (i) is immediate from the spectral decomposition for C° semigroups 
at eigenvalues of finite multiplicity, and the fact that Xj are distinct for |£| > sufficiently 
small, by (H3). Substituting (|3.5p into (|3.ip and computing 

(3.8) 6y(Z, x)=J2 e 2mkx d y {C + 2^) = ^ e 2 ™ kx e^*- 2 ™^ = e^ [8 y {x)i 
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where the second and third equalities follow from the fact that the Fourier transform of 
either the continuous or discrete the delta-function is unity, we obtain 

yielding (|3.7p (ii) by (|3.6p (i) and the fact that <p is supported on [— tt, tt]. □ 

We now state our main result for this section, which uses the spectral representation of 
G 1 and [G|] described in Proposition 13. 21 to decompose the low- frequency Green kernel into 
a leading order piece (corresponding to translational modulation) plus a faster decaying 
residual. Underlying this decomposition is the fundamental relation 

(3.9) G(x,t; y) = £ J e^-^G^, t; yi )}<%, 

which serves as the crux of the low- frequency analysis both here and in [OZ2| IJZ3| . 

Proposition 3.3. Under assumptions (Hl)-(H^) and (Dl)-(DS'), the low-frequency Green 
function G 1 (x, t; y) of (|3.5p decomposes as G 1 = E + G 1 , 

(3.10) E = U'(x)e(x,t;y), 
where, for some C > 0, all t > 0, 

aup||G / (- 1 t,;i/)|| I ,, (!e) <C(l + t)-5 (1 4) 
y 

(3.11) sup W y &{-^ ; y)\\ LP(x) , sup \\d r t G\;t, ; y)\\ LV(x) < C(l + t)'^ 1 '^ 

y y 

sup \\&(; t, ; y)(0, if \\ LP{x) < C(l + t)"^ 1 -?)-* 
y 

forp>2, 1 < r < 2, 

(3.12) S upK^ e (-,^;y)IUn-)^ c '( 1 + *r^ ) " i ^ 1 ^ 

V 

for p>2, < j,l, j + 1 < K + 1, 1 < r < 2, and 

(3.13) sup||^ e (-,t,;y)|| iP(!r) < C(l + 1)-^ 1 "*)-^ 

y 

for < j, I, j + I < K + 1, provided that p > 2 and j + I > 1 or p = oo. Moreover, 
e(x, t; y) = /or t < 1. 
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Remark 3.4. The crucial new observation in the nonconservative case treated here is 
(|3.1ip (iii). which asserts that sources entering in the nonconservative second coordinate 
of the linearized equations experience decay equivalent to that of a differentiated source 
entering in the first coordinate. This is what allows us to treat non-divergence-form source 
terms arising in the second equation of the eventual perturbation equations. 

Proof. Recalling (I3.7P and Lemma [2TTT we have 

i r 2 

G I (x,t;y) = (—) J e^'^UiO^E^QiiMOj^yT^ 
(3-14) " R J= l 



j,k,l=l 



the fact that j3j t i = 0(£ 1 ) suggests the k = 1 terms (corresponding to translation) domi- 
nate the low-frequency Green kernel. With this motivation, we define 

(3.15) e(x,t;y) = (^) J e*^)^) £ M&fafaMZ, y)*dt 



so that 

G J (x, t; y) - U'(x)e(x, t; y) = 

(3.16) ^ u 

+ (^) / ^^^(OE^^ftv^^ 1 ) " U'ix^v^y)*^, 

3,1 

where, by analyticity of v%, vi(£,x) — U'{x) = 0(|£|), and so, by Lemma [2TT1 

(3.17) fofofafox) - U'ix^v^y)* = O(l) 
and 

(3.18) j , 2 p j>l v 2 ^,x)v l ^,yr = O(l). 

Note further that v\ = (1,0) T unless I = 1, in which case /3ji = 0(\£\) by Lemma [2. H hence 

(3.19) dy^jj^v^x) - U'ix^v^y)*) = 0(|e|), 



(3.20) 



(/3 J ,A/(^i(e,x)-c/'(x)) ? )K^yr)(o,i) T = o(|c|), 
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and 
(3.21) 



h[Pj,203MZ> x MZ,vY)=o(\z\) 



(3.22) 



From representation (|3.16p . bounds (|3.17p - ([3.18p . and 3iAj(£) < — 9\£\ 2 , we obtain by 
the triangle inequality 



(3.23) ||G (',t;-)\\L°°(x,y) 



\G l - U'e\ 



to)<C\\e^M\(£)\\»®<CQ. + t)-*. 



Derivative bounds follow similarly, since x-derivatives falling on vjk are harmless, whereas, 
by (|3.19p - (|3.2ip . y- or i-derivatives falling on Vji or on e^'^ x ~ y ^ bring down a factor of |£| 
improving the decay rate by factor (1 + £) -1 / 2 . (Note that |£| is bounded because of the 
cutoff function (ft, so there is no singularity at t = 0.) 

To obtain the corresponding bounds for p = 2, we note that (|3.14p may be viewed 
itself as a Bloch decomposition with respect to variable z := x — y, with y appearing as a 
parameter. Recalling (|l,16p . we may thus estimate 
(3.24) 

sup||G 7 (-,i;y) - U'e(-,t;y)\\ L 2 (x) < 
y 

C SU P HiO^^ki-, zi)vf{-, y)v t {-, y)* \\l^{^{ Zi e[q,x})) 



+ C V] sup 



m^( M ' iX }~ nx) w-,vY 



£ 2 «;£ 2 (*ie[0,X])) 



<C ^2 SU P H0(€) e e|e|2 *llL2(c)Sup||fA : (-,zi)|| i 2 (0iX) ||{; i (-,y)*|| L o 



j,k?i,l 



+ CV S up||0(e)e-^l 2t | 
i,i y 
< C(l + t)-«, 



£2(0 SU P 



lei 



£ 2 (0,X) 



'(0,X) 



Il°°(o,x) 



where we have used in a crucial way the boundedness of v\ in L°° jf| and also the boundedness 
of 

'Vn(^x) - U'(x)\ 

J ~ d(.v n {r) 



9 This is clear for £ = 0, since Vj are linear combinations of genuine and generalized eigenfunctions, which 
are solutions of the homogeneous or inhomogeneous eigenvalue ODE. More generally, note that the resolvent 
of — 7 gains one derivative, hence the total eigenprojection, as a contour integral of the resolvent, does 
too- now, use the one-dimensional Sobolev inequality for periodic boundary conditions to bound the £°° 
difference from the mean by the (bounded) H 1 norm, then bound the mean by the L 1 norm, which is 
controlled by the L 2 norm. 
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in L 2 , where r E (0, £). Derivative bounds follow similarly as above, noting that y- or t- 
derivatives bring down a factor £, while x-derivatives are harmless, to obtain an additional 
factor of (1 + 1) -1 / 2 decay. Finally, bounds for 2 < p < oo follow by L p -interpolation. 
Now, defining 

(3.25) e(x, t; y) := x(*)e(x, t; y), 

where e is defined in (|3. 15|) and x is a smooth cutoff function such that x(i) = 1 for £ > 2 
and x(t) = for t < 1, and setting G := G — U'(x)e(x, t; y), we readily obtain the estimates 
(|3.1ip by combining the above estimates on G 1 — Ue with bound f|3.4j) on G . 

Finally, recalling, by Lemma 12.11 that v\ = constant for I ^ 1 while f5j t \ = 0(|£|), we 
have 

Bounds (13.12P thus follow from (I3.15j) by the argument used to prove (13.11j) , together with 
the observation that x- or t-derivatives bring down factors of ^. Bounds ()3.13j) follow 
similarly for j + I > 1, in which case the integrand on the righthand side of A3. 15H (now 
differentiated in x and or t) is Lebesgue integrable. 

In the critical case j = I = 0, taking t without loss of generality > 1, expanding 

^(o = -^ j -b J e+o(f), 

and setting A(£) := —i^aj — bj£ 2 , we may write e(x,t;y) in (I3.15P as 
(^)j^E4i(0)^ l2 (0)i^(0,y)*e«-^r 1 ^ K) *de 

(3.26) = (^)P-V-jf E 4i (0)^,2 (0)5 2 (0, j^Mr^^df 
= E ^(0)42(0)^(0, y)*(^)P-V.jf e^^r 1 ^^, 

where /3j,i(0) := lim£_>o(£/?j,i(£))> an d the above series is convergent by the alternating 
series test, plus a negligible error term 

(-W. f e l ^ x -^mO{e- e ^ 2t )di 

for which the integrand is Lebesgue integrable, hence, by the previous argument, obeys the 
bounds for j + I = 1. (Note that the integral on the lefthand side of (|3.26p is absolutely 
convergent by £~ 1 (e -mi £* — e~ W2 ^*) ~ \a\ — a,2\t, becoming conditionally convergent only 
when the integrand is split into different eigenmodes.) 

By (D2), we have a,j real and ?Hbj > 0. Moreover, the operator L, since real- valued, 
has spectrum with complex conjugate symmetry, hence bj is real as well. Observing that 
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is an antiderivative in x of the inverse Fourier transform 



(|3.26p is a sum of errorfunctions 

2 



a Gaussian, we find that the principal part 



(3.27) 



x — y — ad 



^ errfn ( - — ^==P^ , < ) «2 (0, y) , 



hence bounded in L°° as claimed, where Oj denote the characteristic speeds of the Whitham 
averaged system and (on further inspection) ^ ■ Cj = 0. This verifies bound f)3. 13j) in the 
final case j = I = 0, completing the proof. □ 



Remark 3.5. See the proof of Proposition 1.5, [OZ2] . for an essentially equivalent estimate 
from the inverse Laplace transform point of view of the critical contribution (|3.26j) . 

3.3 Final linearized bounds 

Corollary 3.1. Under assumptions (H1)-(H4), (Dl)-(DS'), the Green function G(x,t;y) 
of (jl.lip decomposes as G = E + G, 

(3.28) E = U'(x)e(x,t;y), 

where, for some C > 0, allt> 0, 1 < q < 2 < p < oo, < j, k, I, j + 1 < K + 1, 1 < r < 2, 

< C(l + t)-2( 1 /9-l/p)| 



+oo 



(3.29) 



oo 
+oo 



G(x,t;y)f(y)dy 
d r y G(x,t;y)f{y)dy 



LP(x) 



—oo 
+oo 



LP ( X ) 



< C(l + t)-5( 1 /9-l/p)-^ 



d r t G(x,t;y)f(y)dy 



< C(l+t)-2 



l(l/ g _ 1 /p)_| ||r||wnflar+i! 



LP(x) 



+oo 



G(x,t-y)(0,lff(y)dy 



< Cil + t) 



LP(x) 



(3.30) 



+oo 



did*e(x,t;y)f(y)dy 



oo 
+oo 



<(! + *)' 



LP 



& x d*d r y e(x,t;y)f(y)dy 



< (l + t)" 



LP 



did*e(x,t,y)(0,iy f(y)dy 



-i(l/ g -i/ p) -ti+M + i 



i (1/g _l/ p) _(2+M 



<(l+t)-i(l/9-Vp)-^ 



L8, 



LP 



Moreover, e(x,t;y) = /or t < 1. 
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Proof. (Case q = 1). From (|3.1ip and the triangle inequality we obtain 

/ G I (x,t;y)f(y)dy < [ sup \\&(., t; y)\\ LP \f(y)\dy < C(l + ty^M 

and similarly for y- and t-derivative estimates, and products with (0, 1) T , which, together 
with (|3.4p . yield (|3.29p . Bounds (|3.30[) follow similarly by the triangle inequality and (|3.12p - 

(Case q = 2). From (|3.17p - (|3.18p . and analyticity of Vj, Vj, we have boundedness from 
L 2 ([0,X]) — > L 2 ([0, X]) of the projection-type operators 

(3.31) / -> (3jJ Jyl (v n (£, x) - U'(x)) (v h f) 
and 

(3.32) / -> Pj,ikfaMM(*l,f) for k + !. 

uniformly with respect to £, from which we obtain by (|3.16p . (|3.25p . and (|1.16p the bound 



(3.33) 



G I (x,t;y)f(y)dy 



L 2 (x) 



< C\\f\\ L 2 (x) , 



for all t > 0, yielding together with (|3.4p the result (|3.29p for p = 2, r = 1. Similarly, by 
boundedness of Uj, [7' in all L p [0, X], we have 

\ / L^lx) 



e^^khiVk^x)^}) 



L°°(x) 



<Ce-^ 2t \\f(Z,-)\\ LHx) , for kj^l, 



C, 9 > 0, yielding by definitions ([3~T6]) . (ET25]) the bound 

"+0O 



G I (x,t;y)f(y)dy 



< 



(3.34) 



L°°(x) 



(— 

V2vr 



C^)e- e l«l 2 *||/(C,.)|| i2(a;) ^ 

( e) e -^i 2 * 



< C 



L\i,x) 



<C(l + t)-i||/|| iW]) , 



hence giving the result for p = oo, r = 0. The result for r = and general 2 < p < oo then 
follows by LP interpolation between p = 2 and p = oo. Derivative bounds 1 < r < 2 follow 
by similar arguments, using (|3.19p ^ (|3.2ip . as do bounds for products with (0, 1) T . Bounds 
()3.30p follow similarly. 

(Case 1 < q < 2). By Riesz-Thorin interpolation between the cases q = 1 and q = 2, 
we obtain the bounds asserted in the general case 1 < g < 2, 2 < p < oo. □ 

Note the close analogy between the bounds of Corollary 13.11 and those obtained in 
[MaZ3, MaZl] for the viscous or relaxation shock wave case. 
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4 Nonlinear stability 

With the bounds of Corollary (|3.ip . nonlinear stability follows by a combination of the 
argument of [JZ3] IJZ4j and modifications introduced in the shock wave case to treat partial 
parabolicity and potential loss of derivatives in the nonlinear iteration scheme [Zl[ IZ3] . 

4.1 Nonlinear perturbation equations 

Given a solution U(x,t) of (|1.2|) . define the nonlinear perturbation variable 

(4.1) v = U-U = U(x + ip(x,t),t)-U(x), 
where 

(4.2) U(x,t) := U(x + ip(x,t),t) 
and if) : R x R — > R is to be chosen later. 

Lemma 4.1. For v, C7 as in (|4. 1[) . (|4.2j) . and |C/| bounded, 

(4.3) I7 t + /(*/)* - (B(U)U X ) X - g(U) = {d t - L) U'(x)if>(x, t) + P + d x R + d t S, 
where 

(4.4) P= (g(U)-g(U))ij x = {^l) T 0{\v\\^ x \), 

(4.5) R := vip t + B(U)(U X + f.)-^- - (b(U) - B(U)) U x ^ x - B(U)v x ^ x , 

1 + ip x \ J 

and 

(4.6) S := -wj> x . 

Proof. By the definition of U in (|4.2|) we have by a straightforward computation 

U t (x,t) = U x (x + ?l;(x,t),t)ip t (x,t) + Ut(x + ip,t) 
f(U (x, t)) x = df{U (x + V(x, t),t))U x (x + ij>, t) ■ (1 + ^(i, t)) 
C^(z,t) = ^(x + V(x,t),t) • (l + Vfefoi)). 

By t/ 4 + df(U)U x - {B(U)U X ) X - g{U) = 0, it follows that 
(4.7) 

U t + f(U) x - (B(U)U X ) X - g(U) = U x 4> t + df(U)U x ^ x - {B{U)U x ) x i, x - (B(U)U X ^ X ) X 

= U X A - U t iP x + g(U)^ x - (B(U)U X ^ X ) X , 
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where it is understood that derivatives of U appearing on the righthand side are evaluated 
at (x + ip(x,t),t). Moreover, by another direct calculation, using L(U'(x)) = 0, we have 

(8 t - L) U'{x)i> = U x ^ t - U t ip x + df(U)U x ^ x - (B(U)U x ) x ip x - (B(U)U x rp x ) x 
= U x ^ t - Utf x + g(U)ij x - (B(tJ)U x ij x ) x . 

Subtracting, and using the facts that, by differentiation of (U + v)(x,t) = U(x + ip(x,t),t), 

(4.8) U x + Vx = U x {\ + U t + Vt = U t + JJ x ip t , 
so that 

(4.9) U x -U x -v x = -(U x + v x )-^- 1 Ut-U t -v t = -(U x + v x ) '"' 



we obtain 

f/t + /(CO* - (B(U)U X ) X - g{U) = (d t - L)U'{x)i> + v x ^ t - v t ip x 

+ {g{U)-g{U))^ x -(B{U)v x ^ x ) 

(4.10) ( ~ - ib 2 

+ [B(U)(U x + v x " ' 



b(l7) - B(t/)) ^ ' 



X I ; 



yielding flOJ by u x T/>t - = (^t)x - (W'x)t- □ 
Corollary 4.2. T/ie nonlinear residual v defined in (14. ip satisfies 

(4.11) v t -Lv = (d t - L) U'(x)i> -Q x + T + P + R x + d t S, 
where P, R, and S are as in Lemma and Q and T are defined by 

Q : = f(U(x + 1>(x,t),t)) - f(U(x)) - df(U(x))v 

(4.12) - (B(U(x + f/f(x,t),t))U x (x + il>(x,t),t) -B{U{x))U x {x)) 

- (B{U)v x + {dB(U)U x ) v) 

and 

(4.13) T := g(U(x + ^{x, t), t)) - g(U(x)) - dg(U(x))v = (0, l) T 0(\v\ 2 ), 

Proof. Taylor expansion comparing f|4. 3|) and Ut + f(U) x — (B(U)U X ) X — g(U) = 0. □ 
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4.2 Cancellation estimate 

Our strategy in writing (|4.1ip is motivated by the following basic cancellation principle. 
Proposition 4.3 ( |HoZ] ) . For any f(y, s)eL p f) C 2 with f(y, 0) = 0, there holds 

(4.14) / [ G(x,t-s;y)(d s -L y )f(y,s)dyds = f(x,t). 



Proof. Integrating the left hand side by parts, we obtain 
(4.15) 

G(x,0;y)f(y,t)dy - f G(x,t;y)f(y,0)dy + J j \d t - L y )*G(x,t - s;y)f(y, s)dy ds. 
Noting that, by duality, 

(d t - L y )*G(x, t-s;y) = 5(x - y)5(t - s), 

5(-) here denoting the Dirac delta-distribution, we find that the third term on the righthand 
side vanishes in ()4.15j) . while, because G(x, 0;y) = 5(x — y), the first term is simply f(x, t). 
The second term vanishes by f(y,0) = 0. □ 

4.3 Nonlinear damping estimate 

The following technical result is a key ingredient in the nonlinear stability analysis that 
follows. Applying Duhamel's principle to (|4.1ip and using Proposition 14.31 yields 



(4.16) 



f 00 

v(x,t)= I G(x,t;y)v (y)dy 

-oo 

ft />oo 



rt poo 

+ G(x,t-s;y)(-Q y + T + Ry + S s )(y,s)dyds + ip(t)U'{x). 

JO J-oo 



Note that terms Q y and S s involve derivatives of v (respectively second derivative in space 
and first derivative in time) of maximal order, hence to close a nonlinear iteration scheme 
based on (|4.16j) would appear to require delicate maximal regularity estimates rather than 
the straightforward ones that we have obtained. Indeed, estimated using the linearized 
bounds of Corollary 13.11 the righthand side appears to lose several degrees of regularity as 
a function from H K — > I? of v. However, the next proposition, adapted from the methods 
of [MaZ4, Z3], shows that higher-order derivatives are slaved to lower-order ones, hence 
derivatives "lost" at the linearized level may be "regained" at the nonlinear level. This 
effectively separates the issues of decay and regularity, allowing us to close a nonlinear 
iteration without the use of maximal regularity estimates or a more complicated quasilinear 
iteration scheme. 

Proposition 4.1. Let vq € H K (K as in (HI)), and suppose that for < t < T, the H K 

norm of v and the H K+1 norms ofipt(',t) and ip x (-,t) remain bounded by a sufficiently small 
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constant. Moreover, suppose that the Froude number F, viscosity v , and velocity derivative 
u x satisfy the amplitude condition uu x < F . Then, there are constants C,9\ > such 
that, for all <t <T, 

(4.17) \\v(t)f HK <Ce~+ t \\v(0)\\%K+C f (11* + IKVt, ^)\\ 2 h k) (s) ds. 



The proof of this result will be given in Appendix [XJ Here, we briefly outline the main 
ideas. First, notice that by subtracting from the equation (|4.7p for U the equation for U, 
we may write the nonlinear perturbation equation as 

v t + (Av) x - (Bv x ) x -Cv = P- Q{v) x + T(y) + U x i> t - U t i/> x 
(4.18) _ / _ . 

+ g{U)^ x -[B{U)U x ^ x 



where A, B, C are as in (|1.12p . P, Q, and T are as in Corollary 14.21 g and B are as in (jl.lOp . 
and it is understood that derivatives of U appearing on the righthand side are evaluated at 
[x + if)(x,t),t). Using (|4.9p to replace U x and Ut respectively by U x + v x — (U x + v x ) 1 ^ 

and Ut + vt — (U x + v x ) j^ , and moving the resulting vttp x term to the lefthand side of 
(|4.18p . we obtain 



(1 + j> x )vt = (Bv x ) x - (Av) x + Cv + P- Q(v) x + T(v) 
(4 ' 19) + {U x + v x ) i> t + g(U)ip x - (b(U) {U x + v x ) 



Define now the Friedrichs symmetrizer 

'1 



(4.20) S := 



5~ 



2 > 



where 5 2 := — A12 = r 3 (F 1 — vu x ). By (jl.2ip . S is a symmetric positive definite sym- 

f- c -l \ 

metrizer for the hyperbolic part of (|4.19j) in the sense that Y*A = I J is a 

symmetric matrix, where A is as in (|1.12|) . Furthermore, to compensate for the lack of 
total parabolicity of the governing equation, here indicated by the presence of a neutral 
eigenspace of the matrix Si?, we introduce the skew-symmetric Kawashima compensator 

(4.21) K:=rj(^ "A < r, « 1 

and note that, in particular, for n > sufficiently small there exists a constant 9 > such 
that ${KA + EB) > 9. 

Now defining the functional 

K 

S[v) := (v, Ev) + J] [(div, Kd x ~ l v) + (div, Edlv)) 
3=1 
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where here (•,•) denotes the standard L 2 (W l ) inner product, we find by a direct but 
lengthy calculation using Sobolev embedding and interpolation to absorb nonlinear and 
intermediate-derivative terms that 

(4.22) d t £(v) < -fliH^x+C^H^ + UAMf H K (Xit) ) 

for some positive constants C, 6\ > 0, so long as ||C/||jjk remains bounded and the quantities 
i] > 0, ||w||jj-k and \\(tpt, , *Px)\\h k (x) remain sufficiently small. By Cauchy-Schwarz and the 
fact that S is positive definite by (|l,2ip . we have £(v) ~ IMI^k for 7] > sufficiently small 
and hence (|4.22p implies 

d t £(v) < -9 1 £{v) + C (\\v\\ 2 L2 + KM x )f H K {x ,t)) 

from which (I4.17P follows by Gronwall's inequality and, again, the equivalence of £(v) and 
\H 2 H k- 

For more details and a complete proof of the key inequality (I4.22P , see Appendix [A] 

Remark 4.2. The condition (jl.21D gives effectively an upper bound on the allowable 
amplitude of the wave, for fixed Froude number and viscosity. It is not clear that this 
has any connection with behavior. Certainly it is needed for our argument structure, and 
perhaps even for the validity of (|4,17p . which is itself convenient but clearly not necessary 
for stability. The resolution of this issue would be very interesting from the standpoint of 
applications, both in this and related contexts. 

Remark 4.3. The Lagrangian formulation appears essential here in order to obtain quanti- 
tative bounds like (11.211) on the amplitude of the wave. One can carry out damping estimates 
for sufficiently small-amplitude waves in Eulerian coordinates by the argument of [MaZ2 
in the shock wave case; however, the large-amplitude argument of [MaZ4], depending on 
global noncharacteristicity of the wave- corresponding here to nonvanishing of u — s, where 
s is wave speed in Eulerian coordinates- together with bounded variation of U x , appears to 
fail irreparably in the periodic case. As we have shown here, the same argument succeeds in 
Lagrangian coordinates provided that the linearized convection matrix A is symmetrizable 
(the meaning of bound (jl.21fl ). For similar observations regarding the advantages for energy 
estimates of the special structure in Lagrangian coordinates, see [TZlj . 

4.4 Integral representation/^-evolution scheme 

Recalling the Duhamel representation (|4.16|) of the perturbation v along with the decom- 
position G = U'(x)e + G of Corollary 13.11 we find that defining ip implicitly as 



(4.23) 



oo 

ip(x,t) = ~l e(x,t;y)U {y)dy 

-oo 

"t /" + 00 



/ / e(x,t - s;y)(P - Qy + T + Ry + S s )(y,s)dyds, 

JO J-oo 
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where e is defined as in (|3.25p . results in the integral representation 
(4.24) 



/oo 
G(x,t;y)v (y)dy 
-oo 
ft roo 



ft roc 

+ G(x,t-s;y)(P-Q y + T + Ry + S s )(y,s)dyds, 

JO J-oo 



for the nonlinear perturbation v; see \Z1\ IMaZ2] for further details. Furthermore, differen- 
tiating (|4.23p with respect to t, and recalling that e(x,s;y) = for s < 1, 



(4.25) 



/oo 
d{d k x e(x,t;y)U (y)dy 
-oo 

"t f + OO 



/ / di^e(x,t-s;y)(P-Q y + T + R y + S s )(y,s)dyds. 

JO J-oo 



Equations (|4.24j) . (14.25P together form a complete system in the variables (v,dfil>,d k ip), 
< j, k < K + 1, from the solution of which we may afterward recover the shift ip via 
(|4.23j) . From the original differential equation (|4.11|) together with (|4.25p . we readily obtain 
short-time existence and continuity with respect to t of solutions (vjipt^x) £ H by a 
standard contraction-mapping argument based on (I4.17h . (14,23p . and (13.301) . 

4.5 Nonlinear iteration 

Associated with the solution (U,ipt,ipx) of integral system (|4.24p - (|4.25p . define 

(4.26) C{t)-= sup ||(v,V*>^*)||jff*(«)(l + s) 1/4 - 

Lemma 4.4. For all t > for which is finite and sufficiently small, some C > 0, and 
Eq := \\Uo\\ L i nH K sufficiently small, 

(4.27) at)<C(E + ((t) 2 ). 

Proof. By (|4,4p - (|4,6p and (|4.12p ~ (|4,13p and corresponding bounds on the derivatives to- 
gether with definition (|4,26p . 
(4.28) 

\\(P,Q,R,S,T)\\ L i nH 2 < || (v,v x ,^ x )f L2 + 11(^,^,^)11^ < C((tf(l +t)"3, 

so long as \ip x \ < IV'sliT 8 ' — C(^) remains small. Applying Corollary 13.11 with q = 1 to 
representations (|4.24p — (|4.25p . we obtain for any 2 < p < oo 

Wvi-ML^^cii + ty^-^Eo 

(4.29) + C((t) 2 [ (l + t- aj-ia/a-i/p) (t - a )-§ (1 + 

<C( J B + C(i) 2 )(l + i)~ |(1_1/p) 
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and 
(4.30) 

ft 

\\^t^ x K;t)\\ W K +1 , P <C(l + t)- l 2E + CC(t) 2 / (l + t- s )-5d-V P )-i/2 (1 + s) -i ds 

J o 

<c(E + at) 2 )(i + t)-^-yp\ 

yielding in particular that \\{tpti ' 1 Px)\\h k + 1 be arbitrarily small, verifying the hypotheses of 
Proposition Using (fCTTl) and ([4T29]) - (|4\30] ). we thus obtain ||t>(-, i)\\H K {x) — C(Eq + 

C(i) 2 )(l + t)~ 2. Combining this with (|4,3U|) . p = 2, rearranging, and recalling definition 
(14T261) . we obtain flH}. □ 

Proof of Theorem \l.l[ By short-time il^ existence theory, \\(v, ipt, iPx)\\h k ^ s continuous so 
long as it remains small, hence £ remains continuous so long as it remains small. By (|4.4p . 
therefore, it follows by continuous induction that < 2CEq for t > 0, if Eq < 1/4C, 
yielding by (j4T2BT) the result (fT22l for p = 2. Applying (14T29"D ~ (14T3T1 . we obtain (021) for 
2 < p < p* for any < oo, with uniform constant C. Taking p* > 4 and estimating 

II^IU ||0||l2, ||i?||i2, ||S|| L2 , ||T|| L2 (t) < 11(^,^,^)11^ <CE (l + t)~ 1 

in place of the weaker (|4.28p . then applying Corollary 13. II with q = 2 we obtain finally (jl.22p 
for 2 < p < oo, by a computation similar (I4.29|) - (14.30p : we omit the details of this final 
bootstrap argument. Estimate (jl.23p then follows using (|3.30p with q = 1, by 

\\m\\LP < ce (i + 1)& + catf f\i + 1 - sy^-^Hi + s )-?ds 

(4.31) Jo 

<C(l + t)^(Eo + C(t) 2 ), 

together with the fact that U(x, t) — U(x) = v(x — ip,t)+U(x) — U(x — i/;), so that \U(-,t) — U\ 
is controlled by the sum of \v\ and \U(x) — U(x — ~ This yields stability for 

| U — U\ L i nH K\t=o sufficiently small, as described in the final line of the theorem. □ 

A Nonlinear energy estimate 

The goal of this appendix is to prove the inequality (14.22[) . which was the key ingredient 
in the nonlinear energy estimate in Proposition (|4.ip . To this end, we write the nonlinear 
perturbation equation (|4.19p for the variable v = (r, u) T as 

(A.l) (1 + tl) x )vt = (Bv x ) x - (Av) x + Cv+ (U x + v x ) i> t + g{U)ip x + N 

where the function M = ftf(v, Ux^x^t) is defined by 

(A.2) N := P - Q(v) x + T(v) - (b{U) (U x + v x ) J^p) 

10 Note that we have gained a necessary one degree of regularity in if), the regularity of i[> being limited 
only by the regularity of the coefficients of the underlying PDE (|1.2|) . 
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where P, Q, and T are defined as in (|4.4p . (|4.12p . and (|4.13p . respectively. The key to the 
analysis is to carefully keep track of the "hyperbolic" (r) and "parabolic" (u) components 
of v separately. We begin by symmetrizing the linearized convection matrix A of (lA.ip by 
introducing the Friedrichs symmetrizer £ defined in (|4.20p as 



noting in particular that it is symmetric positive definite by the amplitude condition (|1.2ip . 
f-c -1 \ 

The fact that T,A = I ^ J is symmetric then yields hyperbolic properties of the 

solution using straightforward energy estimates, integration by parts, and the Friedrichs 
symmetrizer relation 

St{U,SU x ) = ~(U,S x U) 

valid for all self-adjoint operators S G C nxn and U G C n . Furthermore, for convenience we 
provide here a list of the block-structure of the various matrices arising in the forthcoming 
proofs: notice by definition that 

(A.3) B,B..(l % A„,A„„,= (: »), CO.-? » 





* 



* 



which immediately yields 

(A.4) EB, Y,B X , T. X B, T, XX B, T, X B X 

(A.5) Y1A X , Yj x A, Yj xx A, Yj x A x , YjA xx 

and 

(A.6) SC = SC,= (° I 

We will continually refer back to these dictionaries throughout the proofs in this appendix. 

Remark A.l. The apparently special structure leading to (IA.3P is in fact a special case 
of the more general structure pointed out in shared by the equations of 1-D gas 

dynamics, MHD, and viscoelasticity |TZlllBLeZ| when expressed in Lagrangian coordinates. 

Defining the first-order "Friedrichs bilinear form" as 

Fl[vi,v 2 ] ■= («i,S« 2 ) + {d x vi,Hd x V2) 
our first step in proving Proposition (14. ip is to establish the following lemma. 



See conditions (A1)-(A2) of [TZTI [Zl] . 
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Lemma A. 2. Let v(-,0) £ H 1 and suppose that for < t < T , the H 1 norm of v and the 
H 2 norms ofijj x and ipt remain bounded by a sufficiently small constant. Moreover, suppose 
that the amplitude condition (|1.2ip holds. Then we have the first-order "Friedrichs-type" 
estimate 

\dt' Fl ^ V,V ^ ~ ~ ( Vx > w ^ BVx ) ~ (v xx ,wY,Bv xx ) 

^ A7 ^ +Ci (llwllia + ^KHia+ellTxIlia 

+ C2\\ip x \\ H 2 (IMIhi + ll^ll^) 

+ Y (W^Wm + \\^\\m)+ Fi[v,wAT\ 

valid for all <t <T, for some constants > where w := (1 + ipx)' 1 G L°° . 
Proof. First, notice from (jA.ip and the symmetry of E we have 

(v, T,v) = (v, T<vt) 
2 dt X 1 X ' 

= (v, u>£ {{Bv x ) x - (Av x ) + Cv + U x ip t + v x ip t + g(U)i>x + M)) 
where w := (1 + ipx)~ l - Since T>B is symmetric by (|A.4p then, we have 

(f , toSX-BuaOa;) = - ((ioE x )u, Bv x ) - {v x ,wHBv x ) = - (v, ((wT,) x B) x v) - (v x ,wT,Bv x ) 
and similarly 

(v, wY, (Av) x ) = (v, wT,A x v) - - (v, (wEA) x v) . 
Furthermore, assuming that ||?/>t||#2 remains bounded we clearly have the estimate 

(v, wT, ((U x + v x )if) t + g(U)if) x ) ) = (v, wYTJ x ip. t ) - - (v, (wEip t ) x v) + (y, wY,g(U)i) x ) 

<{\Mh + \\Mh + U*\\h) 

which, by using Cauchy-Schwarz, immediately yields the zeroth-order estimate 

~{v,^v) <-(v x ,wEBv x ) + C(\\v\\l 2 + \\Mh + Hx\\l2) + (v,wm) 

for some positive constant C > 0. 
Continuing, we find that 

7^ {v x ,wEv x ) = (v x ,w x E {{Bv x ) x - (Av) x + Cv+ (U x + v x ) ip t + g{U)^x)) 

+ (v x ,wE {{Bv x ) xx - (Av) xx + (Cv) x + ((U x + v x )tpt + g{U)ip x ) x )) 
+ {v x ,E{wM) x ) 
= :h+h + {v x ^(wN) x ) . 
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(v x ,w x T, {(U x + v x )tp t + g{U)ip x )) < \\w x \\ L °° (\\v\\ 2 h i + \\ipt\\ 2 L 2 + H^llla) 



To estimate I±, notice that (]A.4p immediately yields 

(v x , w x E(Bv x ) x ) = (v x ,w x T,B x v x ) + (v x ,w x T,Bv xx ) < \\w x \\ L °°\\u x \\%i. 
and that, similarly, we have the estimates 

(v x ,w x T,(Av) x ) , (v x ,w x ECv) < \\w x \\ L °°\\v\\ 2 H i, 
by (|A.5|) and (|A.6|) , Finally, noting that for ||V>i||L°° bounded we have 

•VI! • 

we see that together these yield the estimate 

h < {\\v\\ 2 h i + \\u xx \\ 2 L 2 + ||V>t|||2 + HVvllLO • 

To obtain the analogous estimate on I2, first notice that (|A.4|) and the boundedness of 
||io||j^oo, together with Young's inequality, imply 

(v x , wS (Bv x ) xx ) = - ((wE) x v x + wEv xx , (Bv x ) x ) 

= - {v xx ,wEBv xx } - (v xx ,wT,B x v x ) 

- (v x , (wE) x B x v x ) - (v x , (wT,) x Bv xx ) 

1" ' 1 2 1 „||„, ||2 



< - {v xx ,wT,Bv xx ) + C\ I -||ux|li2 + ell^xll^ 

for some constants C\,C% > 0, where e > is a sufficiently small constant to be chosen 
later. Similarly, using (|A.5|) and (|A.6P we find that 

(v x ,wH (Av) xx ) = (v x , wY, {A xx v + A x v x )) - - (v x , (wY,A) x v x ) 

(v x ,wE(Cv) x ) < (\\v\\ 2 L 2 + ^\\u x \\ 2 L 2 + e\\T x \\ 2 L 2 
Finally, noting again that HV'tllL 00 is bounded we find that 

(v x ,wll ((U x + v x )ifj t ) x ) < ((e + \\il>xt\\L°°)\\vx\\ 2 L 2 + ^\\^t\\ 2 H i^ - ^ (v x , {wY 1 ip t ) x v x ) 

<(e+ \\w x \\l°° + ||^ x t|U°°) \\vx\\ 2 L 2 + \\u x \\ 2 L 2 + -\\tf>t\\ 2 Hi 

and ^ 

(v x ,wH (g{U)ip x ) x ) < e\\v x \\ 2 L 2 + -||Vv||ffi- 
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Therefore, by choosing H^tH.?? 2 < e, so that HV^Hl 00 < e by Sobolev embedding, and noting 
that lliAr+H r.oo is bounded, we have 

h < - {v xx ,w'EBv xx ) +C (\\v\\ 2 L2 + ~\\ux\\l2 +e\hx\\ 2 L 2 + ^llV'tHir* + "ll^lln^ 

+ [|u>x||x°° {\\v\\ 2 h i + HWaajllia) 

for some constant C > 0, from which the lemma follows by noting that HwxIIl 00 < HV'xIIh 2 
by Sobolev embedding. □ 

From Lemma IA.2I it follows that if the diffusion ££? were positive definite we would 
immediately have the bound 

\j/M < -B\\v\\ 2 Hl +d (H| 2 + H^ll^ + H^l^x) +Ti[v,wAT] 

by using Sobolev embedding and choosing e > 0, H^Hi? 2 sufficiently small, which, up to 
the contribution of the nonlinear residual terms jV, has the form of the inequality stated 
in Proposition (|4.ip : see the proof of Lemma I A. 3 1 below for details on how this calculation 
would proceed. However, the lack of total parabolicity in the governing equation (jl.2p is 
manifested here in the fact that the matrix EB is not positive definite, rather being only 
positive semi-definite with rank one. In order to compensate for this "degenerate diffusion" , 
we introduce the Kawashima compensator K defined in (|4.21j) as 



K :-- 



-1 

1 



where < ij <C 1 is a small parameter which will be determined later. The fact that the 
hyperbolic effects in (lA.lj) can compensate for this degeneracy in the diffusive term T,B is 
the point of the following lemma. 

Lemma A. 3. Assume the amplitude condition (jl.2ip holds. Then for r] > sufficiently 
small, the matrix Y>B + KA is positive definite and, furthermore, the associated bilinear 
form satisfies the coercivity estimate 

(Z,(zb + ka)0> eiW&Wh+vMiWh) 

for some constant 9 > and all £ = (Cii^2) T G L 2 (1R). 

The proof of Lemma IA.3I is based on simple matrix perturbation argument and is omit- 
ted. Defining now the first-order "Kawashima bilinear form" as 

£ l [vi ,v 2 ] := T\ [vt , v 2 ] + {d x vi , Kv 2 ) 

and noting the special structures 

(A.8) KB,KB.- V § J), KA.- V (l ») , KC = , (j ' 

we have the following refinement of the first-order Friedrichs-type estimate in Lemma IA.2I 
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Lemma A. 4. Under the same hypothesis of Lemma (|A.2j) and for n > sufficiently small, 
we have the first-order "Kawashima-type" estimate 

\j t £ i[ v M ^ -0i (Kllli +v\\tx\\ 2 L 2 ) + -J (Ibllia + ll^ll^i + II^IIhO 

/or some constants 9\,C > 0. 

Proof. Using (1A.8P along with arguments similar to those as in Lemma lA.21 we obtain the 
estimate 

-dt (v x , Kv) = (v x , wK ({Bv x ) x - {Av) x + Cv + U x ip t + + g(U)^ x + N) > 

< CV (jMlh + S\\t x \\ 2 l2 + ^\\u x \\ 2 L2 + ^||n,,||| 2 + ^ (||^||| 2 + ihMio) 
- (v x ,wK Av x ) + {v x ,KwM) 

for some positive constants #i,C > 0, and for any 5 > sufficiently small, where we have 
used that || V* II ~ llV'tlli? 1 can be chosen sufficiently small, say of order 0(6). Since 
w £ L°° then, it follows from Lemma lA.21 that 

-d t £i[v, v]<- (v x , w (EB + KA) v x ) - {v xx wY,Bv xx ) 

+ c i ( f? + 1 ) Ml* + (? + -) H^lli 2 + M + \\ T x\\h + i\\ u xx\\h 



C2\\lpx\\m (iMItfi + Ikxxllia) 

+ C 3 (\ + ]) {\\Mm + \M\m) +Si[v,wJ^\ 



By Lemma lA.31 then, we find that for rj > sufficiently small, say < n < t]q, we have the 
estimate 

- (v x ,w (T,B + KA) v x ) - (v xx wT,Bv xx ) < -9 (\\u x \\ 2 L2 + v\\Tx\\ L 2 + H^llia) 
for some constant 9 > 0. Thus, by fixing 5 and choosing e = e{rf) such that 

9 n9 
(A.9) < 5 = — and < e(rj) = 

2Ci 4 

we find that 

+ C7 1 5) + E )||r a! ||| 8 = -^11^11^. 
By subsequently requiring that the free parameter r/ > satisfy 

(A.10) 0<7/<minj— ,r/ 
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we similarly find that 



rj 

+ T Cl ) II u xx 1 1 £2 5; — o 1 1 1 1 £2 • 



5 

from which it follows by the above requirements on the parameters rj, e(n), and 5 that 
\~m^ v M - ~® (ll^Hi 2 + r ?ll r ^llL2 + H^llla) 

+ Ci (jMli 2 + ~IWIi 2 J + c 2\\^x\\m [\\ v \\m + \\u xx \\ 2 L i) 
+ ^{\\Um + Hx\\ 2 m )+£i[v,wM]. 

Next, the Sobolev inequality [Iffsllfa < 1 1 Sxx 1 1 l 2 1 1 5 1 1 l 2 along with Young's inequality implies 
that ^||ux||| 2 < 2 ll^illxa + jrlMli 2 wn i cn > by now choosing ||V>is||#i sufhciently small 
so that — | + C^UV'xII/fi < 0, completes the proof. □ 

Using similar arguments, we can obtain higher order Kawashima-type estimates by 
defining the fe th -order Kawashima bilinear form as 

k 

£ k [vi,v 2 ] : = ( V ,nv) +J2 {(divuKdt 1 ^) + (a^i,sa> 2 )) 

for each k € N. Indeed, the following estimate can be obtained by simply iterating the 
above argument and using the Sobolev inequality H^Hhj < a ll^c +2 flilL 2 +o~ 1 ||5'||l 2 fo r 
a > sufficiently small. 

Lemma A. 5. Let j E N and v(-,0) 6 and suppose that for < t < T , the H J norm 
of v and the H^ +l norms of ijj x and ipt remain bounded by a sufficiently small constant. 
Moreover, suppose that condition (jl.21j) is satisfied. Then for n > sufficiently small, there 
exist constants 9\, C > such that, for all < t < T, 

Id C 

^ _6l i (iKIIzri +^11^11^-0 + ~2 (IMli 2 + ll^tllL + II^IlL) +£j[v,wN\. 

To complete the proof of Proposition 14.11 it remains to estimate the terms £j[v, wAf] 
corresponding to the nonlinear residual terms in the perturbation equation (jA.ip . In par- 
ticular, our goal is to demonstrate that these terms can be absorbed into the bound already 
computed, in the sense that there exists constants C > and < e <C 1 such that 

(A.ll) £j[v,wM] <e{\\u x \\ 2 HJ + IMl^x) +C(H| 2 + \\iP t \\ 2 Hj + II^IIL) • 

To this end, we notice that from (I4.4h we have the identity 

i 



(g(U) - g{U)) ip x = dg(U + Ov)d0j mp x . 
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Using Sobolev embedding then, we can estimate P in H 1 in a straight forward way, using 
that |M|l°° is assumed to be small (say, at most one). Indeed, using the above integral 
representation for P we immediately obtain 



I^IIl 2 ^ IIW^Hl 2 ~ IMIz^l^slU 00 ^ \\v\ 



and similarly 



\Pt. 



x\\L 2 



< 



d 2 g{U + 0v)d0\ {U x + v x )vi> x 



+ 



L-' 



dg(U + Ov)dO (wl> x ) x 



L 2 



^ IIW^Hl 2 + IKV^Hl 2 \\v\\ L °° + \\{v^ x )x\\l 2 

< IblUall^Hifi + [ba;||i^||V , a;||jyi||w|k ao + \\Vx\\l* W^xWh 1 + \\ v \\l°° \\lpxx 

^ (1Mb 1 + IMIffO ll^llifi 
^ II 11^*11^1 • 

From these estimates together with Cauchy-Schwarz it follows that 

Si[v,wP] < ||u|| L 2||P|| L 2 + ||^|| L 2||P|| L 2 + H^lkalK^P^H^ 
£ IMIffill^llff 1 

where, again, we have used the fact that H^H/fi is small. Since we can control the size 
of ip x in H 1 , it follows that the ||Ua;||i2 term above can be absorbed in the sense that the 
above inequality is of form (lA.lljl . 

Using similar arguments, we can express T(v) as 



T(v) = g(U) - g(U) - dg(U)v 
from which we get the estimates 



d 2 g(U + e.sv)ds ) 9v dJB) v. 



\T{v)\\ L 2 < ||«||x»||«||La £ IMIffi 



and, similarly, 



\\ T (v)x\\ L 2 ^ \\v\\ 2 H i + \\vf H i < \\v\\ H i, 
where we have used the facts that |H|#i is small. As above, these estimates readily yield 

£i[v,wT] = (v,T,wT) + (v x ,KwT) + (v x ,E(wT) x ) 
^ IMIl 2 IM|#i + IMMMI^i 



which again absorbs due to the control over v in H 1 . 

To analyze the remaining terms of A/", consider the term ( B(U)(U X + v x ^ — 



present 



at the end of (|A.2|> . Using the representation B(U) = ( J dB(U + 9v)d6) v + B(U) along 



A NONLINEAR ENERGY ESTIMATE 



34 



with the smallness of v in H 1 and ip x in H 1 , it follows that the associated contribution to 
£i[v,wM] can be absorbed so long as the highest-order Friedrich's term 

v x ,z(w(b(U)(U x + v." 



1 + 



can be shown to absorb. Using integration by parts, we have 
(v x ,e(w(b(U)(U x +v x )-^-] ) ) = -((Zv x ) x ,w(b(U)(U x + v x )-. r ''' 



which absorbs by estimates similar to those previously obtained, using the smallness of v 
in if 1 and ip x in i? 2 . 

To estimate the contribution of the final terms of £i[v,wAf], associated with — Q x , first 
write Q := Q\ — Q2 + Q3 where 

Qi := f(U) - f(U) - df(U)v = (J U d 2 f(U + 6sv)ds^ 9v d9^j v 

Q 2 := B{U)U X - B{U)U X - B(U)v x 
Q 3 := (dB(U)U x ) v 

and notice that the contribution of £\[v, w(Q\ + Qz) x \ absorbs using estimates analogous 
to those obtained above for P and T. To illustrate how to handle the contributions of Q2, 
first notice that by (|4.9p we have 

Q2 = ( [ dB(U + 6v)dJd v) (U x + v x ) - B{U) (l) x + v x ) 

which absorbs as above by the smallness of v in H 1 and ip x in H 2 . 

From the above considerations, then, we immediately have the following lemma. 

Lemma A. 6. Under the same hypothesis of Lemma XA.Hh we have the first- order "Kawashima- 
type" estimate 

j t £i[v, v] < -8, (\\u x f Hl + ||r x .||i 2 ) + C (\\v\\ 2 L2 + U x f m + Utfm) 

valid for some constants 9\,C > 0. 

With this H 1 estimate in hand, the analogous H m estimate follows for any m £ N, 
as in the statement of Lemma IA.51 Finally, using one last time the Sobolev inequality 
1 1 3a; lira ~ \\9xx \\l 2 WdWl 2 1 together with Young's inequality, we have completed the proof of 
the key inequality (j4.22j) . from which the proof of the Proposition 14.11 follows. 
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B High-frequency resolvent bounds 

In this appendix, we carry out the high-frequency resolvent bounds needed for the high- 
frequency solution operator bounds of Section 13.11 To begin, write 

L c = e-* x Le* x = dBd -dA + C, 

where d := {d x + ig). Clearly then, the norm WfW^ := \\df\\ L 2 i[0tX]) + \\f\\ L 2([o,X]) is 
equivalent to the usual norm ||/||ij-i([o,x]) f° r £ e [ — vr ; 7r ] bounded. Further, note that, for 
periodic functions /, g on [0, X], we have the usual integration by parts rule 

(B.l) (f,dg) = (-df,g), 

where (•,•) as above denotes the standard L 2 complex inner product on [0, X}. The main 
result of this appendix is then that for |£| bounded away from zero and sufficiently small 
the resolvent operator (A — L^)~ l is uniformly H 1 — > H 1 bounded for 3?(A) = —r\ < —9 < 
for some constant 9 > 0, which is the content of the following lemma. 

Lemma B.l. Under the derivative condition (|1.2ip . there exist constants C, R > and a 
constant 9 > sufficiently small such that for |A| > R and KA < —9, 

( B - 2 ) \H\m([o,x}) < C||(-k£ - *>)™\\m([o,x])- 

Proof. For £, K , as defined in the proof of Proposition ^. 1\ define the first-order Kawashima- 
Bloch bilinear form as 

£i[vi,V2] := {vi,Hvi) + (dvx^Kvz) + (dv\,Y,dv^) 

and suppose w is a solution of (A — L^) w = f. Then using the coercivity estimate of Lemma 
IA.3( it follows by taking the real part of the equation 

£i[w, (A - L{)w] = £i[w, f] 

and using the equivalence of £i[w,w] ~ ||Hlj?i> we obtain similarly as in the proof of 
Proposition 14.11 

(B.3) (KA + 9)\\w\\ 2 Al + e\\B%w\\h < C(\H\h + \\f\\%), * > 0- 

Similarly, defining the analogous first-order Friedrich's-Bloch bilinear form 

Fi[vx,V2] := (ui,Si;i) + (dvi,Hdv^j 

and taking the imaginary part of the equation 

F\ [w, (A - L{) w] = T\ [w, f] 
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we obtain 

(B.4) |3A||H| 2 < C(|Ml|i + \\Bd 2 x w\\ 2 L2 + „, M//J 

Summing (|B.3|) with a sufficiently small multiple of (|B.4p . we obtain for KA > —9/2 

lAiHIi^cdHI^ + ll/lllO, 

yielding the result for |A| > 2C by equivalence of -ff 1 and H . □ 

C The subcharacteristic condition and Hopf bifurcation 

At equilibrium values u = r -1 / 2 > 0, the inviscid version 

i 



(C.l) C7 t + /([/), = ( 7(c/) ) ) g(C/) = l-r^, /(IT) 



of (|1.8p has hyperbolic characteristics equal to the eigenvalues ±^7^ of df, and equilibrium 

characteristic ^- equal to d T f(r, u*(t)), where u*(t) := defined by q{T,u^{r)) = 0. 

The subcharacteristic condition, i.e., the condition that the equilibrium characteristic speed 
lie between the hyperbolic characteristic speeds, is therefore 



u 3 u 3 



(a2 > T < ^ 

or F < 4 as stated in Remark 11.21 

For 2x2 relaxation systems such as the above, the subcharacteristic condition is ex- 
actly the condition that constant solutions be linearly stable, as may be readily verified by 
computing the dispersion relation using the Fourier transform. For the full system (jl.8p 
with viscosity v > 0, a similar computation, Taylor expanding the dispersion relation about 
£ = 0, reveals that constant solutions are stable with respect to low-frequency perturbations 
if and only if the subcharacteristic condition F < 4 is satisfied. 

Next, let us examine the profile ODE cV + ((2F)" 1 r" 2 ) / = 1 - r(q - ct) 2 - cv(t- 2 t')' 

— 1/2 

near an equlibrium uq = (q — ctq) = t q > 0, and examine the circumstances for which 
Hopf bifurcation occurs. Linearizing about r = To, and rearranging, we obtain 

-2_// 1 ( „2 „2\J 1 / n o/2 — c \ „ ._ n 



(C.3) cur, V + (<? - ci)r> + p^— : r, ... 

V Uo/Z / y/F 

for which the eigenvalues are the roots \i of ctfj, 2 + /3^ + 7 = 0, where a = cvt~ 2 , (3 = c 2 — c 2 , 
and 7 = /9 . Considering this as a problem indexed by parameters uo, c, and g, we see 
that Hopf bifurcation occurs when roots /J,j(uo, c, q) cross the imaginary axis as a conjugate 
pair, i.e., when j3 = and 7 > 0. 
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These translate, using (|C.2p to the Hopf bifurcation conditions 



(C.4) c = c s = -^= and F > 4. 

s/F 



Experiments of [Nl| indicate that bifurcation occurs at minimum wave speed, i.e., as c 
increases through the value c s . That is, the minimum wave speed among nontrivial periodic 
waves is 

and the minimum value of F for which nontrivial periodic waves occur is F > 4. The 



frequency at bifurcation is lj = y/jja = t^ 2 u 1 ^ 2 \J {^/F — 2), and the period is X = 
So prescribing X as we do, we must choose F > 4, then solve oj = ^ to obtain 

(C.6) r = ,W ^ 



■X 2 (x/F-2) 

Nearby this value and with c near c s , we should find small-amplitude periodic waves. 

Remark C.l. The above discussion shows in passing that, similarly as observed in the 
conservative case in |OZl| . small-amplitude periodic waves arising through Hopf bifurcation 
from constant solutions are necessarily unstable as solutions of the time-evolutionary PDE, 
since they inherit (a small perturbation of) the necessarily unstable dispersion relation of 
the limiting constant solution from which they bifurcate. On the other hand, in the large- 
amplitude limit, roll waves might well be stable. As observed by Gardner (see [G jlOZlp . this 
is determined by stability of the bounding homoclinic wave, which in the conservative case 
was known to be unstable. A good starting point for the study of roll waves, therefore, might 
be to determine linearized stability of solitary pulse solutions corresponding to homoclinic 
solutions of the profile ODE. Evidence for linearized stability of some viscous roll waves is 
given in [Nl], namely, the approximate Dressier waves arising in the small viscosity limit. 



D Numerical stability investigation 

We conclude by suggesting a number of practical techniques for the numerical testing of 
stability. These can be carried out either in the Eulerian coordinates of [N2] or in the 
Lagrangian coordinates of this paper. As suggested in a more general setting in [B] , several 
of the algorithms may be eaily adapted from an existing nonlinear evolution code. Com- 
parison of these different methods, and determination of stability in different regimes, are 
interesting problems that we hope to carry out in future work. 
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D.l Method one: the power method 

An easy numerical method to approximate the function := max Kcr(L^) determining 
stability, with R(£) < for £ 7^ corresponding to (Dl) and R(£) < — corresponding to 
(D2). (Condition (D3) can be verified by an Evans analysis, as was already done in some 
cases in [N2] and elsewhere.) 

The method is just to approximate numerically the time-evolution of the linearized equa- 
tion wt = L^w on [0,X] with periodic boundary conditions, which should be a straight- 
forward adaptation/simplification of the nonlinear code you have already written to study 
nonlinear stability with respect to periodic perturbations, and which should work quite well. 
Denote the solution operator as e L ^. Then, a good approximation is 



where / is a square wave pulse centered at x = X/2 and T is large, say T = 10, T = 50 
or T = 100. This should be relatively straightforward, and plotting agains £ for 

£ G [—it, 7r] should quickly tell stability. See [BMSZ] for related investigations and discussion. 

D.2 Method two: discretization 

Instead of Evans computations as in |OZl| (these involved finding the zero-level-set of a 
two parameter Evans function, with reported problematic results) one could alternatively 
proceed from a Bloch decomposition/matrix linear algebra point of view. 

That is, one could discretize L^ on [0, X] with periodic boundary conditions as a large 
tri-diagonal matrix 



acting on vectors (Ui, ... ,Ul) of sample points, where Uj ~ U(Xj/L), and virtual point 
Uq = Ul (periodicity), and A is a discrete derivative, for example the forward difference over 
h := X/L, treating [0,X] as a torus to wrap generate needed values Uj for j < or j > L. 
For each £, call the fast linear algebra functions in MATLAB to generate the real part of 
the largest real part eigenvalue of T as a function R(£). If R(£) < and R(£) < — c£ 2 , 
c > 0, then we have spectral stability- otherwise not. This should be fast even for 100 x 100 
matrix or so. The discretization in £ is over [— it, it], so also no problem. 50 points should 
suffice. 

Note [B] that discretizaton of the linearized operator L is typically already done for a 
standard method-of-lines realization of the linearized time evolution. 

Remark D.l. In an interesting recent talk by Dwight Barkley [B] . he pointed out that 
doing power method for e Lt , t small @ with implicit scheme is something like using inverse 
power law on (/ — Lt)^ 1 . Note that (J — Lt)^ 1 is expected to be compact for t small in 
parabolic problems, so this is a preconditioning step paralleling the Fredholm theory or 
Birman-Schwinger approach on the analytical side. 



(D.l) 



RiO^T 1 log 




(D.2) 



T(0 := (A + i£)B(A + iO - (A + i£)A + C, 



This is essentially equivalent to the method suggested in Appendix lD.il 
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D.3 Method three: nonlinear evolution 

The simplest test of course is just to run the full nonlinear problem on a large domain 
[—NX, NX], N >> 1, with periodic boundary conditions and square pulse wave initial 
conditions centered at x = 0. If the difference between the solution and the unperturbed 
periodic wave remains bounded in L°°, then the wave is stable, otherwise not. The exper- 
iment should be run only up to time T « NX to avoid interactions with the boundary. 
This, and the sensitivity of numerical evolution of nonlinear equations, are the main disad- 
vantages of the method. The advantage is that this can be converted from existing nonlinear 
code for evolution on a single period [0, X] (easy to change). A variation is to solve the 
linearized equations vt = Lv := (d x Bd x — d x A + C)v numerically, which would be more 
stable but require modification (straightforward, however) of the nonlinear code, changing 
over to linear. 



D.4 Method four: Evans function computations 

A final approach is to compute the Evans function D(£, A) (straightforward [OZ1} |N2] ) and 
plot zero level sets of -D(£, •) for varying £ (harder). This is not recommended in the basic 
form just described- in practice this was time-intensive and gave poorly resolved results 
|OZlj . A somewhat more reasonable variation would be to plot just the level sets near 
(£, A) = (0,0) (difficult, due to crossing/singularity at the origin, but contained) to verify 
(D2), then use winding number computations for D(£, •) to verify (Dl). 
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